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Abstract 

We construct a new Hartree-Fock-Bogoliubov (HFB) mass model, labeled HFB-18, with a 
generalized Skyrme force. The additional terms that we have introduced into the force are 
density-dependent generalizations of the usual t\ and t% terms, and are chosen in such a way as 
to avoid the high-density ferromagnetic instability of neutron stars that is a general feature of 
conventional Skyrme forces, and in particular of the Skyrme forces underlying all the HFB mass 
models that we have developed in the past. The remaining parameters of the model are then 
fitted to essentially all the available mass data, an rms deviation a of 0.585 MeV being obtained. 
The new model thus gives almost as good a mass fit as our best-fit model HFB-17 (a = 0.581 
MeV), and has the advantage of avoiding the ferromagnetic collapse of neutron stars. 

PACS numbers: 21.10.Dr, 21.30.-x, 21.60.Jz, 26.60.Dd, 26.60.Kp 
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I. INTRODUCTION 



Astrophysical considerations require that one have available nuclear-mass models as rig- 
orously based as possible. In this way one might hope to be able to extrapolate from the 
mass data, which cluster fairly closely to the stability line, out towards the neutron drip 
line, and make reliable estimates of the masses of nuclei so neutron rich that there is no 
hope of measuring them in the foreseeable future; such nuclei are found in the outer crusts 
of neutron stars, and also play a vital role in the r-process of nucleosynthesis. To this end 
we have developed a series of nuclear-mass models based on the Hartree-Fock-Bogoliubov 
(HFB) method with Skyrme and contact-pairing forces, together with phenomenological 
Wigner terms and correction terms for the spurious collective energy; all the model param- 
eters are fitted to essentially all the experimental mass data (see Ref. l| and references 
quoted therein). To make the extrapolations to neutron-rich nuclei as reliable as possible, 
the underlying Skyrme force in model HFB-9 [2] and all later models was constrained to fit 
the zero-temperature equation of state (EoS) of neutron matter (NeuM), as calculated by 
Friedman and Pandharipande [3] (FP) for realistic two- and three-nucleon forces. (We have 
so far been unable to obtain mass fits as good as our published ones when constraining to 
the more complete, and slightly stiffer, realistic neutron-matter EoS labeled A18 + 5v + 
UIX* [4]. Our findings are consistent with a recent analysis of the n~ /n + ratio in central 
heavy- ion collisions indicating that this EoS is too stiff [5|.) 



Because of the neutron-matter constraint our models can be used to extrapolate beyond 
the neutron drip line and calculate with some confidence the EoS of the inner crust of neutron 
stars |6j (throughout this paper we assume zero temperature). Since the good agreement 
of our forces with the FP calculation [3J of neutron matter extends to the highest density 
encountered in neutron stars it might be thought that our inner-crust EoS is continuous with 
that of the homogeneous core, the transition between the two regions taking place at around 
0.5po ; where po (— 0.16 nucleons.fm -3 ) is the equilibrium density of symmetric nuclear 
matter (SNM). However, in fitting our forces to the neutron matter of FP [3], we assume 
that our ground state is spin unpolarized, but in fact the underlying Skyrme forces of all our 
models, like all conventional Skyrme forces of the form (pQ), predict that beyond a certain 
supernuclear density the ground state of NeuM becomes ferromagnetic, i.e., at least partially 
polarized [3, [si, [q]. In the case of the Skyrme force BSkl7, the force underlying our best-fit 
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model, HFB-17 [K, complete polarization sets in at a density of Pf rm g — 1.24 po- On the 



other hand, microscopic calculations using different realistic forces and different methods 



11 



12 



13 



10 



14] all predict no such polarization, at least up to about 5p - Moreover, in the case 



of all our own previously published HFB forces the predicted ferromagnetic state is unstable 
against collapse, the energy becoming more and more negative as the density increases 
(see, for example, the lowest curve in Fig. [fl constructed for force BSkl7). This predicted 
ferromagnetic collapse sets in at densities that are certainly encountered in the cores of all 
neutron stars, and is contradicted by the very existence of neutron stars. 

Actually, the core of neutron stars does not consist of pure NeuM but rather of so-called 
neutron-star matter (N*M), which just below the crust is composed of neutrons with a weak 
admixture of proton-electron pairs in beta equilibrium (for simplicity we will assume that 
these are the only particles present also at higher densities). But if NeuM were indeed un- 
stable against collapse then N*M would be likewise, since, at the very least, it would always 
be energetically advantageous for N*M to spontaneously transform into NeuM through elec- 
tron capture beyond some critical density. Thus the stability of NeuM against collapse is 
a necessary condition for the existence of neutron stars. But even if no such instability is 
implied, there is a general tendency for Skyrme forces of the conventional form (pQ) to lead 
to the ground state of N*M being polarized. However, calculations based on realistic forces 



15|- 



show the ground state of N*M, like that of NeuM, to be unpolarized at all densities 

Our purpose here is to show that by adding suitable terms to the Skyrme force it is 
possible to eliminate the anomalous prediction of a ferromagnetic transition in neutron stars, 
with essentially no impact on the high-quality fits to the mass data that we have previously 



obtained with conventional Skyrme 
been followed by Margueron et al. 



orces. (An alternative approach to this problem has 



16 



17|.) In Section II we discuss in more detail the 



nature of the spurious transition to a ferromagnetic state in NeuM and N*M associated 
with conventional Skyrme forces of the form (CQ), considering not only our own BSkl7 [l|, 
but also the widely used SLy4 [lsj], which was specifically constructed for neutron-star 
calculations. Section III shows how extra terms in the Skyrme force can stop this spurious 
transition, while in Section IV we describe the new mass fit. Our conclusions are summarized 
in Section V, and in the Appendix we present the full formalism for the generalized form of 
Skyrme force used here. 
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II. FERROMAGNETIC INSTABILITY 



The Skyrme forces that we have used in all our previously published HFB models have 
the conventional form 

v itj = t (l + x P CT )5(ry) + + xiP^)^ [p| Sfa) + <5(r -) pj] 
+t 2 (l + x 2 P a )^p ij .8(r ij )p ij + ~t 3 (l + x 3 P (J )p(r) a 5(r Ji ) 
+^2^0(0'* + o-j) .p^ x 8(rij)Pij , (1) 

where r^- = — r^, r = (r^ +rj)/2, p^ = —ih(Vi — Vj)/2 is the relative momentum, P CT 
is the two-body spin-exchange operator, and p(r) = p n (r) + p p (r) is the total local density, 
p n (r) and p p (r) being the neutron and proton densities, respectively. 

The fit of the parameters of this form of force to the nuclear masses, and the various 
other constraints mentioned above, leave us absolutely no freedom to avoid an unphysical 
ferromagnetic collapse of NeuM. The situation for our force BSkl7 is shown in Figs. [Hand [2j 



(The corresponding curves for force SLy4 18| shown in these two figures are discussed below.) 
These curves have been calculated using Eq. (C. 14) of Bender et al. which gives the 
energy per nucleon of homogeneous nuclear matter with arbitrary charge asymmetry and 
degree of polarization; we define the latter quantity for nucleons of charge type q (q = n or 

p) b y 

h q = P -^± • (2) 

The curve labeled "BSkl7: no polarization" in Fig.fUshows the energy per neutron calculated 
under the constraint I an = 0, while the curve "BSkl7: polarization allowed" shows the same 
quantity calculated at each density by minimizing the energy per neutron with respect to 
I an . The first of these curves is seen to be in excellent agreement with the realistic EoS of 
NeuM given by FP [3], but the second curve shows that this agreement is destroyed when 
polarization is allowed, the system collapsing. The corresponding value of I an for BSkl7 is 
shown in Fig. [2j the rapidity with which polarization sets in with increasing density will be 
seen. 

As for N*M, to the nuclear energy corresponding to the Skyrme force in this neutron- 
proton mixture we have to add the electron kinetic energy, for the density of which we take 
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the exact expression (see, for example, Section 24.6c of Ref. [201 ]) 



24tt 2 V h ) 
where 



{-8x 3 + 3x(l + 2x 2 )(l + x 2 ) 1/2 -3sinh" 1 x} , (3) 



x — 

mc 



— {^ Pe ) ll \ (4) 



p e being the electron density (equal to the proton density p p , because of electric-charge 
neutrality). We minimize the total energy per nucleon with respect to the proton fraction 
Ye = Pp/ Pi while imposing the constraint I an = I ap = 0, and obtain the curve labeled 
"BSkl7: no polarization" in Fig. [31 If we next minimize the total energy with respect to 
Y e , I an and I ap we obtain the curve labeled "BSkl7: polarization allowed" in Fig. [31 the 
corresponding values of Y e , I an and function of density are shown in Fig. HI It will 

be seen that once polarization is admitted BSkl7 leads in N*M to the same instability with 
respect to collapse that we found in NeuM (Fig. [1]). 

The case of Skyrme force SLy4- Figs. [T] and [3] also show the corresponding curves for 
force SLy4 ]18j in NeuM and N*M, respectively. It will be seen that in NeuM (Fig. [[]), 
the onset of polarization has been postponed to pjrmg = 4.4p , and that although there is 
a considerable softening of the EoS there is no collapse, at least at neutron-star densities. 
However, Fig. [3] shows that although there is still no collapse in N*M, the transition to a 
spurious ferromagnetic state, and the associated softening of the EoS, takes place at the 
much lower density of Pfrmg = 2.5po, which certainly lies within the range of densities found 
in neutron stars. This shows that stability against a ferromagnetic transition in NeuM at a 
given density does not guarantee stability in N*M at the same density. (A similar conclusion 
has been reached in Ref. [li]].) 

The fact that SLy4 has greater stability against polarization than does BSkl7 is a result 
of setting x<i = — 1. Now within the framework of a conventional Skyrme force of the form ([1]) 
we have been unable to find an acceptable mass fit for X2 = — 1, and indeed SLy4 performs 
badly as a global mass model, the rms deviation for the even-even nuclei being quoted as 



5.1 MeV [2l(]. This could have implications for the composition of both the outer and inner 
crusts, and might explain the differences between SLy4 and BSkl4 predictions shown in 
Figures 1 and 2 of Ref. [6J]. It seems that it is impossible to have both a good mass fit and 
stability against polarization with the conventional form (Tj[|) of Skyrme force. 



III. STABILITY RESTORED 



We note now that the form ([I]) does not exhaust the possibilities for a Skyrme-type force, 
and we shall consider here two extra terms, writing our complete Skyrme force as 

<i = v l , 3 + ^t i (l + x 4 P a )^{ P y(r)i 3 5(r l] ) + 5(r ij )p(rf P 2 lJ } 

+ t 5 (l + x 5 P a )^ PlJ .p(ry5(r t3 ) PlJ . (5) 

The £4 and t$ terms are density-dependent generalizations of the t\ and £ 2 terms, respectively. 
The formalism for this generalized Skyrme force is developed in Appendix |A], where we 



show in particular that Eq. (C. 14) of Bender et al. [19j for homogeneous nuclear matter of 
arbitrary charge asymmetry and polarization can be generalized to include the new terms 
simply by making the substitutions of Eqs. flA31al) - flA31dl) . 



We consider now just the simpler case of NeuM, the stabilization of which is a necessary 
condition. For complete polarization in the presence of the conventional Skyrme force 
the energy per nucleon is given by 



3/r , 



mirn (6vr 2 p) 2/3 + ^(6vr 2 ) 2 / 3 t 2 (l + x 2 )p 5 / 3 . (6) 

The only Skyrme term operative here is the one in i 2 . The ferromagnetic collapse of neutron 
matter predicted by all our Skyrme forces at supernuclear densities arises from the fact that 
the combination t 2 (l + ^2) is always negative. However, stability can always be enforced by 
adding a new repulsive term in £5, as can be seen from Eqs. f)A31cj) and flA31dl) . But we 
do not want to disturb in any way the unpolarized configuration of NeuM, since we know 
from the experience with our recent mass models that with the conventional form of Skyrme 
force ([1]) alone it is easy to fit the realistic EoS of FP jjj, and wish to do so here with the 
new model. As the energy per nucleon of this latter configuration in the presence of the 
conventional Skyrme force ([1]) is given by 

e = I ^-(37r 2 p) 2/3 + ^o(l-xo)p + ^3(l-x 3 K +1 

+ A( 37 r 2 ) 2 / 3 {£ 1 (l-x 1 ) + 3t 2 (l + x 2 )}p 5 / 3 (7) 

(see Eq. (IA26f) ). it follows from Eqs. (1 A3 lap and ()A31b[) that the t$ term can be completely 
canceled in unpolarized NeuM by adding a £4 term with its parameters constrained by 

= 7 (8) 



and 



U(l-x A ) = -3t 5 (l + x 5 ) . (9) 

It is now highly convenient to require that the stabilizing terms in t± and £5 cancel com- 
pletely in unpolarized nuclear matter of any charge asymmetry. This leads, using Eq. (1A13j) . 
to a second condition on t 4 , 

£ 4 = -3*5(5 + 4x5) , (10) 

which, combined with Eq. leads to 

4 + 5x 5 

x 4 = — . (11) 

5 + 4x 5 K J 

Thus all three parameters of the £4 term will be completely determined by the parameters of 
the £5 term, i.e., £5,0:5 and 7. These latter three parameters leave us with ample flexibility 
for stabilizing NeuM against polarization, and possible collapse. We stress, however, that it 
will also be necessary to check the stability of N*M. 

IV. THE HFB-18 MASS MODEL 

Even though the new terms exactly cancel in homogeneous nuclear matter they will not 
do so in finite nuclei, and we cannot simply add them on to the BSkl7 force. Rather it will 
be necessary to make a complete refit of all the model parameters to the mass data. 

Our HFB calculations for finite nuclei are performed exactly as for the HFB-17 model lj. 
In particular, the treatment of pairing, which is neglected in the neutron-matter constraints 
discussed in the previous section, is highly realistic. As usual, we take a contact pairing 
force that acts only between nucleons of the same charge state q, 

vf r ( ri , rj ) = v^[p n (r),p p (r)] 5{r l3 ) , (12) 

where the strength v nq [p n , p p ] is a functional of both the neutron and proton densities. But 
instead of postulating a simple functional form for the density dependence, as is usually 
done, we construct the pairing force by solving the HFB equations in uniform asymmet- 
ric nuclear matter with the appropriate neutron and proton densities, requiring that the 
resulting gap reproduce exactly, as a function of density, the microscopic 1 5'o pairing gap 



calculated with realistic forces 



221 ]. We follow our usual practice of allowing the proton pair- 



ing strength to be different from the neutron pairing strength, and for allowing each of these 
strengths to depend on whether there is an even or odd number of nucleons of the charge 
type in question. These extra degrees of freedom are taken into account by multiplying 
the value of v nq [p n , p p ], as determined by the nuclear- matter calculations that we have just 
described, with renormalizing factors f^, where /+, /~ and /~ are free, density-independent 
parameters to be included in the mass fit, and we set = 1. 

To the HFB energy calculated for the Skyrme and pairing forces we add a Wigner cor- 
rection, 

E w = V w ex V {-\(^^\ \+V^\N-Z\expl- (^\ J , (13) 

which contributes significantly only for light nuclei with N close to Z. Our treatment of 
this correction is purely phenomenological, although physical interpretations of each of the 



two terms can be made 



23. 



A second correction that must be made is to subtract from the HFB energy an estimate 
for the spurious collective energy. As described in Ref. 24J ]. the form we adopt here is 

E cM = EZt nk {b tanh(c|/3 2 |) + d\{3 2 \ exp{-Z(|/3 2 | - /3 2 ) 2 }} , (14) 



in which E^ nk denotes the cranking-model value of the rotational correction and /3 2 the 
quadrupole deformation, while all other parameters are free fitting parameters. 

The final correction that we make is to drop Coulomb exchange. This is a device that 
we have successfully adopted in our most recent models, beginning with HFB-15 [25], and 
it can be interpreted as simulating neglected effects such as Coulomb correlations, charge- 
symmetry breaking of the nuclear forces, and vacuum polarization. 

The parameters of the model, i.e., of the Skyrme and pairing forces, and of the Wigner 
and collective corrections, are fitted to the same set of mass data as are all our models since 
HFB-9 |2j, namely, the 2149 measured masses of nuclei with N and Z > 8 given in the 2003 
Atomic Mass Evaluation [2y]. This fit is subject to the constraints on both unpolarized 
and polarized NeuM discussed in the previous section, as well as our usual requirement that 
the isoscalar effective mass M* take the realistic value of 0.8M in SNM at the equilibrium 
density po- The values of the Skyrme, pairing and Wigner parameters resulting from this 
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fit are shown in Table fl] (ea is the pairing cutoff parameter 

HQ) 

: this defines the BSkl8 



"force" . The parameters of the collective correction are s 



The starting point for the new fit was the force BSkl7 [l], and we stress that the parameter 
search was far from exhaustive, particularly with regards to £5, 25 and 7. Our choice of these 
latter parameters was somewhat arbitrary, being limited only by the requirements that the 
energy-density curve for polarized NeuM lie entirely above that for unpolarized NeuM, and 
that there be no significant deterioration in the quality of the overall mass fit. That the 
first of these requirements is satisfied is seen in Fig. [5j which automatically guarantees the 
stability of NeuM against polarization. Proceeding as in Section II, we have checked also 
that N*M is stable against polarization over the same density range, as seen in Figs. Eland [7] 
(force BSkl7 gives virtually identical results in N*M, provided we do not allow the ground 
state to be polarized for this force). Likewise, it is apparent in the first line of Table II III 
that we have satisfied the requirement on the quality of the mass fit, and indeed it can be 
seen from this table that in some respects the new model is better than the old. The fact 
that large regions of the extended parameter space remain unexplored allows ample scope 
for future improvement, both with respect to the quality of the mass fit and conformity of 
the force to reality. Note, however, that as long as our force is constrained by Eqs. ([8]), ffTUl) 
and (fTTj) the £4 and £5 components of the force will not contribute to the effective mass, and 
it will not be possible to exploit the advantages of a surface-peaked effective mass {27]: see 
Eq. flATOl) . 

The comparison of the new parameter set with the original BSkl7 set in Table [I] shows 
that the changes in the t\ and £2 terms are much greater than those in the to and £3 terms; 
this simply reflects the fact that the new t\ and £5 terms cancel in homogeneous unpolarized 
matter, and act only through the gradient terms to which they give rise. (Their contribution 
to the binding is nevertheless far from negligible, amounting to 38 MeV in the case of 208 Pb.) 
For the same reason there is very little difference between the BSkl7 and BSkl8 values of 
the droplet-model parameters |28j shown in the first five lines of Table IIV1 There is likewise 
very little change in most of the Landau parameters shown in the last eight lines of Table HVl 
Only Gq and G\ show any substantial differences between BSkl7 and BSkl8, these being 
the only two Landau parameters to which the £4 and £5 components can contribute when the 
condition ([8]) holds and the particular choice £4 = t 5 is made. In particular, the generalized 
Skyrme force BSkl8 leads to Go being larger than -1 in SNM for all densities, thus removing 



lown in Table HT1 
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the spin instability that was predicted by our previous force BSkl7, as can be seen in 
Figure [9j On the other hand, a spin-isospin instability still occurs in SNM at about the 
same high density as that found for BSkl7, the two forces giving very similar values of G' 
for all densities (see Figure ITUT) . However, insofar as beta-equilibrated N*M has, as we have 
seen, a large neutron excess, this spin-isospin instability in SNM is of no consequence for 
neutron stars. 

Finally, we have checked that with the new interaction, BSkl8, causality is satisfied in 
NeuM for all densities found in neutron stars, i.e, the velocity of sound v s is smaller than 
the velocity of light c. The former is given by the relativistic expression 



*-*/§3S • (15) 



in which P(p) is the pressure, 



ode(p) 



P{P) = P 2 -^ (16) 
op 



and £{p) is the total energy density, 

£(p)=p(e(p)+M n c 2 ) , (17) 



e(p) being simply the energy per neutron, as given by Eq. flA26j) . The calculated value of 
v s /c for BSkl8 is shown as a function of p in Fig. [HI 

Using this interaction BSkl8, we have constructed a complete mass table, labeled HFB- 
18, running from one drip line to the other over the range Z and N > 8 and Z < 110. The 
results are very similar to those for HFB-17, the rms difference between all 8389 predictions 
being 0.433 MeV, and the mean difference (HBF-18 - HFB-17) 0.198 MeV. We have also 
calculated the spins and parities of all these nuclei; only for about 10% of these nuclei is 
there a difference with respect to the HFB-17 prediction. 



V. CONCLUSIONS. 



We have extended our earlier Skyrme-HFB mass models by the inclusion of terms that 
are density-dependent generalizations of the usual t\ and £2 terms. We have shown that 
these new terms can be chosen in such a way as to prevent the high-density ferromagnetic 
collapse of neutron stars that was a general feature of our previous HFB mass models, 
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without compromising the excellent fit to the mass data that we obtained in the past, and 
without relaxing any of the previously imposed constraints of conformity to reality. The 
mass predictions made by the new model, HFB-18, are, in fact, very similar to those made 
by the preceding model, HFB-17 1|. These two models not only give better fits to the mass 
data than does any other published model except that of Duflo and Zucker 29j, but they 
are also by far the most microscopically founded models, and in particular their underlying 
interactions (BSkl7 and BSkl8, respectively) have been fitted to realistic calculations of 
both the EoS and the 1 S'o gap of neutron matter. They can thus be expected to make more 
reliable predictions of the highly neutron-rich nuclei that appear in the outer crust of neutron 
stars and that are involved in the r-process. Moreover, these mass models can be used to 
extrapolate beyond the drip line to the inner crust of neutron stars, using the respective 
interactions to calculate the EoS in this region. Our confidence in this extrapolation derives 
not only from the fit of the interactions to neutron matter but also from the precision fit to 
masses, which means that the presence of protons and the existence of inhomogeneities in 
the inner crust are well represented. 

Finally, pursuing our aim of a unified treatment of the different regions of a neutron star, 
we can use these effective interactions to calculate the EoS of the homogeneous core (at least 
in its outer parts where no complications from the possible appearance of hyperons and other 
particles arise). Of course, for pure NeuM nothing new can be obtained in this way beyond 
what has already been given by the realistic calculations to which our effective interactions 
have been fitted, but these interactions can then be reliably used to calculate N*M, which 
is not treated in the realistic calculations of FP 3j. Another important application of these 
effective interactions, which would hardly be practical with realistic forces, is to make a 
detailed study of the transition between the inner crust and the fluid core. 

We have seen that for N*M, as for masses, the two interactions BSkl7 and BSkl8 give 
very similar results, provided we assume that the ground state for the former is unpolarized, 
which in fact is not the case. Here lies the basic difference between the two interactions: with 
BSkl8 the ground state of NeuM and N*M is indeed unpolarized at all densities prevailing 
in neutron stars. This elimination of the spurious polarization is the essential contribution 
of this paper, made possible by the introduction of the terms in £ 4 and t 5 . However, we 
have so far made only a partial search in the new, extended parameter space, and there 
remains the possibility not only of improving the fit to the mass data still further but also 
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of imposing further realistic constraints on the Skyrme force. 



APPENDIX A: FORMALISM FOR GENERALIZED SKYRME FORCE. 

The formalism for the conventional Skyrme force ([I]), with expressions for the energy den- 
sity, single-particle fields, etc., has been given many times, and is conveniently summarized 
in Brack et al. jsoj ]. The extension to cover the t A terms was given by Farine et al. [27], but 
the t$ terms have not, to our knowledge, been dealt with before, except for the special case 



3l| 



where the exponent 7 is set equal to 1 
Energy density. 

Assuming invariance under time reversal, the HFB energy is written as the integral of a 
purely local energy-density functional 



#hfb = / ^HFB(r) d 3 r 



(Al) 



where 



£hfb(t) = ^Sky 

+ £coul 



p n (r), Vp„(r), r n (r), 3 n (r),p p {r), Vp p (r), r p (r), J p (r) 



Pp{r) 



£- 



pair 



Pn {r),p n {r),p p {r),p p {r) 



(A2) 
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The first term here, the energy density for the Skyrme force of this paper, is given by 



S 



Sky 



= £ 

q=n,p 
1 



1 

+ 4 h 



2M n Tq + 2 to 



1 + ) P z 



1 + 



Q+*o) E^ 

^ij (pr + ^Vp) 2 ) - Q + *i) E + 



q=n,p 



1 + pr - -(Vp) 



^(vp) 2 ) + Q+x 2 ) e (p^-^vp,) 2 ) 

' q=n,p ^ ' 



+ ^ 



1 + I P J 



1 

+ ^4 



+ 17*4 



1 

+ I' 5 



x 7 q=n,p 

1 + ^ (pr + Jj(Vp) 2 ) - Q + z 4 ) E \Pi r i + i( V ^ 



q=n,p 



1 + ^4 J p(Vp) 2 - Q + X 4 J Vp • E P<Fpq 

' \ / q =n p 



1 + 2 X s ) ( P T 



q=n,p 



J 2 

<-l 



- ^ {tixi + t 2 x 2 ) J 2 + ^ (h - t 2 ) E J < 

q=n,p 

- ^ (t 4 x 4 p /3 + hx 5 p<) J 2 + Yq ( Up& - hpl ) E J < 
+ \w I J-Vp+ E -VVpJ . 



q=n,p 



(A3) 



It will be seen that the £5 term has the same form as the t 2 term, multiplied by the p 7 
factor. No such simple relation between the t 4 and t\ terms is apparent, because we have 
eliminated terms containing a Laplacian, through integration by parts over the entire system 
(this accounts for the terms linear in (3). 

The second term in Eq. (1A2I) is the Coulomb energy density, which, since we are neglecting 
Coulomb exchange (see Section IV), is given simply by 

1 



£coul — 2" e Pch^ 



Coul 



(A4) 



in which ep cn is the charge density associated with protons (this differs from ep p because 



we are taking account of the finite size of the proton 



24|), and V Conl is the electrostatic 
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potential, given by 



V 



Coul/ 



e / drr 



, Pch{r') 
\r — r' 



(A5) 



The last term in Eq. (1A2I) is the pairing-energy density, discussed fully in Refs. Q . 

Self- consistent s.p. fields. In coordinate-space (assuming time-reversal invariance), the 
HFB equations read 



E 



ti q (r) aa > - \ q 8 a „> 



<t'=±1 



-h' q {r) a(Tl + \8 aal 



Ei 



A q (r)5 aa , 

where the s.p. Hamiltonian h' (r) aa i and pairing field A g (r) are given by 



(A6) 



h 2 

h' q {r)a'a = -V • 2M , (r) VcW + U q {r)8 aa , - iW q (r) ■ V x a a , a 



(A7) 



and 



The s.p. fields appearing in Eq. ( 1A7I) are defined by 



(A8) 



h 2 



2M*(r) 



U q {r) 



dr q (r) 



d Pq (r) 



*p£L W,(r) = (A9) 



From Eq. ( 1A3I) . we find 



2M* 



h -ti 

2M ? 4 



1 



1 



1 

4 2 



1 



1 + -x 2 ) P + - + x 2 ) p q 



1 



1 

+ -U 



1 + -x 4 I p - I - + x 4 ) p q 



1 + o X 5 P + ^ + X 5 ) Pq 



p\ 



(A10) 
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U q = t 



1 

+ 4 f > 



1 

+ -h 
4 



1+ 2 X °) p ~ \\ + X ° ) /) " 



1 + 5*0 (-+^V) + Q+x 2 ) (r, + iW, 



+ — *3 

12 6 



1 

+ n*4 



(l + l -x^j (2 + a)p« +1 - Q + a*) ^2p% + ap*" 1 £ p 2 , 

1 + la*) p^ 1 {2(l + /3)pr - (2/9 + 3) (^(Vp) 2 + pV 2 p) } 
+ U + ff^fopVp • Vp 9 + 3p 2 V 2 p g - 2p 2 r ? 

+ /3(/3 - l)p 9 (Vp) 2 + /3pp 9 V 2 p - ^/3p [( V ^') 2 + 4 ^'V - 2 /vV%'] } 

+ ^5 J {(1 + l)pr + ^7(Vp) 2 + \p^ 2 p) 
+ (^ + ^5) {pr 9 + ipV 2 p 9 + 7 ^ \p q >T q '-~(Vp gl ) 2 \ + -jVp-V Pg } 

\ J " q'=n,p ^ > 

- ^(t 4 ^A 1 + WV 2 + ^MA 1 -t57P 7 - 1 ) E J J 



q—n,p 



4 ^ «9p g 

q'=n,p * 



(All) 



and 



W 9 = Jw V(p + p g ) --(t 1 x 1 + t 2 x 2 )J+-(t 1 -t 2 )J q 

Z o o 

- -^{Ux^pP + t 5 x 5 p~ / )J + -(t 4 / - ^P 7 )^ • 



(A12) 



The Coulomb field V Cou \r) is given by Eq. (jA"5l) . 

Unpolarized homogeneous nuclear matter. For the energy per nucleon of infinite nuclear 
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matter of density p and asymmetry r\ (defined by 77 = (p n — p p )/p), Eq. (1A3I) reduces to 

-^{^ (1 + " )5/3 + ^ (1 -" ,5/3 } 

3 - (2x„ + 



+ \to 



3 

+ 40 (1 



3 

+ — 1 2 
40 



+ s t3 



3 

+ 40< 4 



3 

H *5 

40 



(2 + xi) J P 6/3 (77)-(- + xi)F 8/3 (77) 



(2 + x 2 )F 5/ 3(r ? ) + (- + x 2 )F 8/3 (r ? ) 



P&F 



3-(l + 2x 3 )r7 2 



(2 + x 4 )F 5/3 (r / )-(- + x 4 )F 8/3 (r ? ) 



(2 + x 5 )F 5/3 (r ? ) + (- + x 5 )F 8/3 (r ? ) 



n~< +1 k 2 
p tv F , 



(A13) 



where 



(A14) 



and 



(l + ri)'+(l-r,)' 



(A15) 



For SNM (77 = 0) Eq. flAl"3|) reduces to 



e(r) = 0) 



3ti 2 



k 2 F + -t p + ^77 



10M * 8 
9 



80 



3ti + t 2 (5 + 4x 2 



P4 + ^3P Q+1 



+ + — t 5 (5 + 4x 5 )p 7+1 A; 



SO 



80 



F 5 



(A16) 



where 



2 1 1 

M ~ JU + W v 



(A17) 



16 



From Eq. (1A17j) we have 
B 2 e 

k° 



F dkl 



rj=0 



3ff 
5M 

3 



9 

k 2 F + -t p + 



3ti + t 2 (5 + 4x 2 ) 



pk\ 



+ 




16 




9 


+ 




80 




3 


+ 




80 



(A18) 

1 /3 

Evaluating this at p = po an d introducing kpo = (37r 2 p /2) defines the usual incompress- 
ibility 



d 2 e 



hpo dkl 



9p; 



rj=0, p=po 



2 d*e 



(A19) 



rj=0, p=p 



the two forms being equivalent only because at p = po we have 

= . 



de 
dp 



(A20) 



v=o, p=Po 

Expanding (1 ± rj) x up to p 2 allows us to rewrite Eq. (1A13|) as 

;2 



e(v = 0) + \k 2 F (±.-^ v + e symV 2 + O (rf 



(A21) 



where we have introduced the symmetry energy 



j sym 



(P) 



1L 

6M 



k 2 F --t (2x + l)p + 



1 

24 



-3tixi + t 2 (4 + 5x 2 ) 



pk 2 F 



-^t 3 (l + 2x 3 )p a+1 - l -t 4 x 4 pP +l k 2 F + ^t 5 (4 + 5x 5 )p^ +1 k% 



(A22) 



(note the breaking of charge symmetry implied by the neutron-proton mass difference; this 
is explicitly included in the finite-nucleus calculations). The usual symmetry coefficient 28] 
is then given by 



J 



(po 



(A23) 



For the density-symmetry coefficient 28] defined by 

de 



L 



de 

Kr 



dk f 



3p - 



sym 



P=P0 



dp 



(A24) 



P=P0 
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it follows from Eq. (1A22|) that 



h 2 



+ 



3M 

a + 1 

16 

5 + 3 7 

24 



k F0 ~ o*o (2x + l)p + 7T7 



24 



t 3 (l + 2x 3 K +1 -^— ^ 2 



t 5 (4 + 5x 5 )p2 +1 4o 



-3tiXi + t 2 (4 + 5x 2 ) 



Po k F0 



(A25) 



Setting i] = 1 in Eq. (1A13j) gives us for the energy per nucleon in unpolarized NeuM 
3h 2 13 9 

+ lt s (l - x 3 )p Q+1 + lt 4 (l - ^)P" +1 k Fn + ^ 5 (1 + * 5 )p 7+1 4„ , (A26) 



where 



k 



Fn 



(3vr 2 p) 



2 nW 3 



(A27) 



(Eq. flA26l) cannot be derived from Eqs. flMB and flA22l l 

Note that the terms in t 4 and t 5 will not contribute to any of the foregoing expressions 
for homogeneous unpolarized matter if the conditions (IE]), ( ITOl) and ( TTTT) are satisfied. 

Polarized nuclear matter 

In general the Skyrme energy density can be decomposed into a time-even part gi ven 
by Eq. (1A3|) and a time-odd part which is non-zero only when time-reversal invariance 
is not satisfied. In polarized homogeneous nuclear matter these two parts take the form 

/r r q + Cgp 2 + C{ ( Pn - p v f + Clpr + C{{p n - p p ){r n - r p ) , (A28) 



c even 
C Sky 



and 



<?odd 
C Sky 



2M r 



CqS 2 + Cj(s n Sp) 2 + CqS ■ T + (s n Sp) • (T n 



T I 



(A29) 



where s 9 and T q are the spin density and kinetic spin density respectively (for a precise 
definition, see, for example, Bender et al. |l9]), and s = s n + s p , T = T n + T p . 

The expressions for the coefficients appearing in Eqs. (IA28P and (IA29I) in the case of the 
conventional Skyrme force (CE]) can be found, for example, in Appendix B of Ref. 19]. The 
coefficients Cq, Cf , Cq and Cf depend only on the t and t 3 terms of the Skyrme force ([I]) 
and therefore remain unchanged when the new terms of Eq. (T5]) are included. This is not the 
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case for the other coefficients. However Cq and C\ can be readily obtained by comparing 
Eqs. (IA28j) and flA3j) . The expressions for the remaining coefficients Cq and CT can also be 
obtained from Eq. f]A3|) using the gauge invariance of the Skyrme force 0, 34]: — Cq and 



— Cj coincide with the coefficients in front of the terms proportional to J 2 and (J n — J p 



(A30a) 



in Eq. f1A3h . The various coefficients are thus given by 

= -t + — t 3 p a 



(A30b) 



(A30c) 



C{ 



(A30d) 



CI 



3 1/5 \ 3 R 1 /5 
16 1 4 2 U 2 / 16 H ^4 5 \4 51 ' 



(A30e) 



(A30g) 



(A30h) 



Under the constraints (jBj), (Ej) and ( 1101) . the new terms in the Skyrme force do not affect 
the coefficients Cq and G\. In the particular case of model HFB-18 for which £4 = t&, the 
coefficient Cf is also unchanged, so that only Cq changes. 
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It can thus be seen that the expressions of the above C-coefficients given in Ref. 19j can be 
extended to the generalized Skyrme force (jSJ) simply by making the following substitutions: 



U — > U + Up 



P 



(A31a) 



t\X\ — > t\X\ + t^x^pr 



(A31b) 



t 2 ^t 2 + Up 1 , 



(A31c) 



t 2 X 2 — > ^2 + ^5^5P° 



(A31d) 



With this simple rule, the expression for the energy per nucleon of asymmetric polarized 
homogeneous nuclear matter can be easily obtained from Eq. (C.14) of Ref. 
Landau parameters 

The dimensionless Landau parameters for symmetric nuclear matter corresponding to 
the generalized Skyrme force (j5J) are given in terms of the C-coefficients of Eqs. ( lA30al) - 
( lA30hj) bv 



F = N 



2C£ + 2G5k* + 4p^° + p 2 ^ + pr^f + ^°(2r + 2pkf 



dp dp 2 



dp 2 dp 



(A32a) 



F^ = N 



*2C^ ~~\~ QiC^ /up 



F\ = -2N Q C T kl 



F[ = -2N Cl4 



G = N 



2Cq -\- 2Cq k-p 



(A32b) 

(A32c) 
(A32d) 

(A32e) 



G' = N 



(A32f) 



d = -2N C^k 2 F 



(A32g) 
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G[ = -2N Q C{kl (A32h) 

where Nq is the density of s.p. states at the Fermi level 

2M*k F , K x 

n it 2 

For the conventional Skyrme force (CD), the above expressions reduce to those given in Ap- 



pendix D of Ref . [19( . Since with BSkl8 Cq is the only C-coefficient that is modified by the 
t 4 and t 5 terms it follows that they affect only the Landau parameters G and G\. 
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FIG. 1: (Color online) Energy per nucleon (T= 0) for pure neutron matter (NeuM) with forces 
BSkl7 and SLy4. The squares are the results of the realistic calculations from Friedman and 
Pandharipande Q]. 
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FIG. 2: Neutron polarization I an for pure neutron matter (NeuM) with force BSkl7. 
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TABLE I: Force BSkl8: lines 1-16 show the Skyrme parameters, lines 17-21 the pairing parameters, 
and the last four lines the Wigner parameters (see text for further details). For convenience of 
comparison we also show model HFB-17 [1] (note that in the latter reference there were misprints 
in the pairing parameters f^ 1 )- 





HFB-18 


HFB-17 


to [McV fm 3 ] 


-1837.96 


-1837.33 


t\ [MeV fm 5 ] 


428.880 


389.102 


ti [McV fm 5 ] 


-3.23704 


-3.17417 


t 3 [MeV fm 3+3a ] 


11528.9 


11523.8 


t 4 [MeV fm 5 + 3 ' 3 ] 


-400.000 


- 


t 5 [MeV fm 5+3 T] 


-400.000 


- 


X 


0.421290 


0.411377 


Xx 


-0.907175 


-0.832102 


X2 


57.7185 


49.4875 


X3 


0.683926 


0.654962 


X4 


-2.00000 


- 


%5 


-2.00000 


- 


Wq [MeV fm 5 ] 


138.904 


145.885 


a 


0.3 


0.3 


P 


1.0 


- 


7 


1.0 


- 


f + 


1.00 


1.00 


fn 


1.06 


1.05 


U + 


1.04 


1.04 


fp 


1.09 


1.10 


EA [MeV] 


16.0 


16.0 


Vw [MeV] 


-2.10 


-2.00 


A 


340 


320 


Vyf [MeV] 


0.74 


0.86 


Ao 


28 


28 
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TABLE II: Parameters of Eq. (|14p for collective correction to model HFB-18. 



b (MeV) 


0.8 


c 


10.0 


d (MeV) 


3.0 


I 


16.0 


P°2 


0.1 



TABLE III: Rms (<r) and mean (e) deviations between data and predictions for model HFB-18; for 
convenience of comparison we also show model HFB-17 [1]. The first pair of lines refers to all the 



2149 measured masses M that were fitted 



261 ]. the second pair to the masses M nr of the subset of 



185 neutron-rich nuclei with S n < 5.0 MeV, the third pair to the neutron separation energies S n 
(1988 measured values), the fourth pair to beta-decay energies Qp (1868 measured values) and the 
fifth pair to charge radii (782 measured values 32|). The last line shows the calculated neutron-skin 
thickness of 208 Pb for these models. 





HFB-18 


HFB-17 


<j(M) [MeV] 


0.585 


0.581 


e(M) [MeV] 


0.007 


-0.019 


a(M nr ) [MeV] 


0.758 


0.729 


e(M nr ) [MeV] 


0.172 


0.119 


a(S n ) [MeV] 


0.487 


0.506 


e(S n ) [MeV] 


-0.012 


-0.010 


<r(Qp) I Me V] 


0.561 


0.583 


e{Qp) [MeV] 


0.025 


0.022 


cr(R c ) [fm] 


0.0274 


0.0300 


e(Rc) [fm] 


0.0016 


-0.0114 


^( 208 Pb) [fm] 


0.15 


0.15 
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TABLE IV: Parameters of infinite nuclear matter for force BSkl8; for convenience of comparison 
we also show force BSkl7 





BSkl8 


BSkl7 


a v [MeV] 


-16.063 


-16.054 


PO [fm" 3 ] 


0.1586 


0.1586 


J [MeV] 


30.0 


30.0 


K v [MeV] 


241.8 


241.7 


L [MeV] 


36.21 


36.28 


M*/M 


0.80 


0.80 


M*/M 


0.79 


0.78 


F 


-0.12 


-0.12 


K 


0.97 


0.97 


Fi 


-0.60 


-0.60 


Fl' 


0.032 


0.068 


Go 


-0.33 


-0.69 


G 


0.46 


0.50 


Gi 


1.23 


1.55 


G\ 


0.50 


0.45 
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— BSkl7: no polarization 
BSkl7: polarization allowed 

— SLy4: no polarization 
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FIG. 3: (Color online) Energy per nucleon (T= 0) for neutron-star matter (N*M) with forces 
BSkl7 and SLy4. 
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FIG. 4: (Color online) Proton fraction Y e , neutron and proton polarizations l an and I^p, respec- 
tively, for neutron-star matter (N*M) with force BSkl7. 
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FIG. 5: (Color online) Energy per nucleon (T=0) for pure neutron matter (NeuM), polarized 
and unpolarized, with force BSkl8. The squares show the results of the realistic calculations of 
Friedman and Pandharipande [3j]. 



200 I i | i | i | ' | ' | ' | ' | i | i | r 




FIG. 6: Energy per nucleon (T= 0) for neutron-star matter (N*M) with force BSkl8. System 
stable against polarization in ground state. 
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FIG. 7: Proton fraction Y e for neutron-star matter (N*M) with force BSkl8. 
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FIG. 8: Speed of sound v s (in units of the light speed c) for pure neutron-matter (NeuM) with 
force BSkl8. 



29 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

P [fm 3 ] 

FIG. 9: Landau parameter Go in symmetric nuclear matter for forces BSkl7 (dashed line) and 
BSkl8 (solid line). 
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FIG. 10: Landau parameter G' Q in symmetric nuclear matter for forces BSkl7 (dashed line) and 
BSkl8 (solid line). 



30 



